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We present an intuitive formalism for implementing cellular automata on arbitrary topologies. 
By that means, we identify a symmetry operation in the class of elementary cellular automata. 
Moreover, we determine the subset of topologically sensitive elementary cellular automata and find 
that the overall number of complex patterns decreases under increasing neighborhood size in regular 
graphs. As exemplary applications, we apply the formalism to complex networks and compare the 
potential of scale-free graphs and metabolic networks to generate complex dynamics. 



Introduction — Cellular automata (CA) on graphs in 
principle provide the possibility to monitor systematic 
changes of dynamics under variation of network topol- 
ogy. In practice, however, unambiguously studying 
the relation between topology and dynamics with CA 
is conceptually difficult, since changes in topology in- 
evitably induce changes in the rule space. Proposed by 
von Neumann [1 as a model system for biological self- 
reproduction, a surge of research activity from the 80's 
onwards [2 established them as the standard tool of com- 
plex systems theory and spatio-temporal pattern forma- 
tion [3 on regular grids. Another discrete and binary 
modeling approach for complex biological systems are 
random Boolean networks (RBNs), introduced by Kauff- 
man [4]. While the CA framework introduces one rule 
for all regularly ordered cells with bi-directional links, 
the original RBNs consist of randomly and directionally 
linked nodes with individual rules. Here, we present a for- 
malism that generalizes CA to arbitrary architectures. It 
allows (i) the establishment of a general correspondence 
between CA and isotropic RBNs and (ii) the comparison 
of the potential of different topologies to generate com- 
plex dynamics. As applications we examine the topolog- 
ical sensitivity of elementary CA, monitor the number of 
complex rules of CA under increasing neighborhood size, 
and compare the dynamic potential of scale-free graphs 
and representations of metabolism as substrate graphs. 

The formalism — Within the CA framework, the discrete 
(binary) state Xi G S = {0, 1} of a node i at time 
t -\- 1 solely depends on its own state and the states 
of its d neighboring nodes at time t. All cells are up- 
dated synchronously by the same, time-independent rule 
/ : ^d+i _^ ^ implement CA on a directed or undi- 
rected graph G, we have to account for different neigh- 
borhood sizes di due to the heterogeneous connectivity 
and thus, in general, to allow for individual rules 
Our strategy instead is to impose constraints on the rule 
space, motivated by simple physical requirements, in or- 
der to obtain a set of discrete rules, implementable on 



arbitrary topologies: 

• Homogeneity = / V z, i.e. the same rule applies 
to all nodes in the graph. 

• Isotropy / = f{xi^ pi)^ i.e. rules may not depend on 
the order of neighboring states and are thus func- 
tions of the density of neighboring states, pi{t) = 
i" ^ij^ji^)' Here, G is represented by the adja- 
cency matrix A: If a link connects node j to node 
i, Aij = 1, and we call j an input node of i. The 
number of all input nodes is called the in-degree of 
node i, di = Aij. 

• Functional simplicity, i.e. the rule / is a piecewise 
constant function of the density pi. 

Elementary Cellular Automata — The simplest CA, 
termed elementary CA (EC A) [2 , are defined on a one- 
dimensional grid with minimal neighborhood size, d = 2, 
and a binary state space, S = {0, 1}. The 2^ = 8 dif- 
ferent neighborhood configurations result 
in 2^ = 256 possible rules. In this set, 2^ = 64 rules 
fulfill the conditions mentioned above and depend only 
on the state Xi (0 or 1) and on the density pi of neigh- 
boring states (0, 1/2, or 1). These 64 rules are called 
outer-totalistic [2 and are now parametrized with the 
rule parameter set (a, /3, 7): 
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We distinguish the following cases for the rule parameters 
7: The state Xi{t + 1) may be or 1 independently 
of the state Xi{t) itself, or it may remain unchanged (+) 
or be flipped (— ), a,/3,7 G {0,1,+,—}. The frequently 
used majority rule [3 El |7l H] , for example, where a node 
i is mapped onto 0/1 if the density pi is below/above 
0.5, and stays in its state otherwise, is described in our 
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formalism by (a,/3, 7) = (0, +,1). For a,/3,7 G {0,1}, 
the corresponding C A rules are called totalistic [2] , since 
Xi{t-\-l) depends exclusively on the density pi of the input 
states. Only these rules have strict RBN rule equivalents 
(see Table |!|. 

Aside from the initial system state x(0) := 
(xi, X2, . . . , xat) at t = 0, the patterns of rule (0,0,0) 
and rule (1,1,1) are perfectly symmetric under the 
action of the operator T : ^ 1 — ^ G {0, 1}. 
The operator T exchanges all Os and Is in an array of 
elements, which can be both a pattern consisting of 
O's and I's or a set of rule parameters. Note that the 
elements {+, — } remain unaffected under the action of 
T. Generally, the symmetric rule to (a, 7) is rule 
T(7,/3, a). The patterns emerging from the action of 
a rule onto an initial state, written as (a,/3, 7) • x(0), 
are identical to the inverted patterns emerging from 
the inverted initial state Tx(0) due to T(7,/3, a): 
(a,/3,7) • x(0) = r(7,/3,a) • Tx(0). Explicitly, the 
symmetric rule to (0,1,+), corresponding to the EGA 
with rule number 218 |2], is (+,0,1) with EGA rule 
number 164 (see Table [l| for more examples). Some rules, 
like the majority rule (0, +, 1), are self-symmetric. After 
elimination of all symmetric counterparts, 34 different 
EGA rules remain. 
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FIG. 1: Rules (1,0,-), (+,+,-), (1,+,-), and (1,0,+) (in 
EGA terms 37, 108, 109, and 133), on graphs with N = 100 
nodes and d = 2. (a) The upper row shows patterns from 
a regular EGA architecture, the lower row from randomized 
counterparts. Time runs downwards, (b) The Shannon en- 
tropy S changes differently for every one of the four rules 
(colored lines) under randomization and is not trivially corre- 
lated with the variation of the diameter of the graph (dotted 
black line), as a prominent topological observable. The results 
are qualitatively independent of the network size. 



= 1 (P? log2 + pj log2 pj ) . 

^ denote the ratios of O's 



Which of these 34 rules are topologically sensitive? That 
is, which lead to patterns of considerably different com- 
plexity when implemented on the regular EGA grid and 
the RBN architecture? Wolfram classified GA heuristi- 
cally according to the complexity of the emerging pat- 
terns into the four Wolfram classes [9 . On graphs, this 
classification inevitably fails because of a lacking natural 
node order. Instead, we apply two entropy-like measures, 
the Shannon entropy S and the word entropy which 
we have previously shown to provide a feasible framework 
for the quantification of pattern complexity [lOl [TTJ [12] . 
The Shannon entropy S serves as a measure for the ho- 
mogeneity of the spatio-temporal pattern, by averaging 
over all nodes: S = Y^iL 
The probabilities p^ and p, 
and I's in the time series of node i. The word en- 
tropy W serves as a complexity measure beyond sin- 
gle time steps. It quantifies the irregularity of a time 
series by counting the number of words, i.e. blocks of 
constant states confined by the respective different state: 

W = ^ Ez^i (- ELi^ilog2Pi)- The probability p\ is 
the number of words of length / divided by the number of 
all words found in the time series of node i. The maximal 
possible word length is given by the length t of the time 
series analyzed. 

We compare 10^ random initial conditions on the regu- 
lar architecture with 10^ samples of a randomized graph, 
where the number of incoming and outgoing links of ev- 
ery node is preserved and kept to (i = 2, but the link 
architecture has been randomized [13 . Notably, for a 
considerably large number of randomization steps, we 
generate random regular graphs [14] rather than Poisson- 
distributed random graphs [T5^. 

We classify rules as topologically sensitive, if the differ- 
ence of the mean entropies for regular and randomized 
architectures is beyond the standard deviation of the dif- 
ference. Out of the 34 rules, 14 rules fulfill that condition 
for at least one observable, S or W. These can be divided 
into three groups: (1) For some rules, the ratio of con- 
stant and oscillating nodes changes under randomization, 
but complex or chaotic patterns never occurs. (2) Oth- 
ers exhibit chaotic patterns on both topologies, but the 
amount of complexity varies. (3) The most interesting 
rules are those, for which the change of topology leads 
to a fundamental change in the complexity of the result- 
ing patterns, i.e., a change in the Wolfram class. These 
rules are (1,0,-), (+,+,-), (1,+,-), and (1,0,+) (in 
EGA terms rules 37, 108, 109, and 133). Typical time 
evolutions of these four rules on regular and random ar- 
chitectures are shown in Figure [ija). Moreover, these 
four rules react specifically to topological changes. Fig- 
ure [TJb) shows the Shannon entropy against the number 
of randomization steps performed. While the patterns 
of rules (1,0,—) and (+,+,—) change already when a 
small number of shortcuts are introduced into the system. 
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(1, +, — ) stays constant in this regime but shows higher 
S and large variations for strongly disordered topologies. 
Finally, the Shannon entropy of the patterns emerging 
from rule (1,0,+) grows monotonously with the random- 
ization depth. 

Regular graphs — How much complexity is possible on reg- 
ular graphs? With growing neighborhood size (i, the 
number of possible densities p and therefore the num- 
ber of possible rules increases. For the sake of simplicity, 
we restrict our investigation to a binary state space and 
to rules with a single threshold: k,: 



a , pi < n 
(3 , Pi = 
7 , pi> K 



(2) 



For networks with d = 12, 11 different threshold param- 
eters K e {1/12, 2/12,..., 11/12} lead to 336 different 
rules, where symmetric rules are considered only once. 
To estimate the number of rules with complex (that is 
in our context: non- trivial) patterns, we calculate W for 
time evolutions. The word entropy is a feasible complex- 
ity measure for individual time evolutions. It however 
fails to disentangle periodic patterns (Wolfram class II) 
from complex (Wolfram class IV) ones. We therefore 
supplement our classification with a detrended fluctua- 
tion analysis (DFA) [16 . This method characterizes the 
time correlations of a signal with a single scaling expo- 
nent, by calculating the variance of the signal from its 
trend in a time window for different window sizes. The 
DFA exponent is the slope of the mean variance against 
the window size and lies between 0.5 and 1.5 for white 
and Brownian noise, respectively. Applied to the time 
evolution of the system's state density p{t) — ^iXi(t) 
as, e.g., in [7], it can be used to discriminate stationary 
and periodic patterns from complex ones. We count pat- 
terns as complex, if the DFA exponent is positive and 
the word entropy > 1. However, the exact values of 



these thresholds do not alter our results qualitatively. As 
shown in Figure [2j the number of possible rules (dashed 
line) increases linearly with while a maximum of com- 
plex patterns (full lines) occurs for d = A. The striking 
overall reduction of complexity for neighborhood enlarge- 
ment, seen as the most dominant effect in Figure |2] can 
be understood qualitatively from a homogeneity ratio- 
nale: In the limit of a fully connected graph, all nodes 
see nearly the same neighborhood and thus follow the 
same dynamics (see ^] for a more detailed explanation 
of a corresponding phenomenon). 

Complex networks — The formalism of Eq. ([2| can be 
transferred to networks of arbitrary topology. Compared 
to regular graphs with global neighborhood size the 
case where pi = n will rarely occur in graphs with hetero- 
geneous connectivity. We thus simplify the set of possible 
rules with single threshold: By setting a = (3 m Eq. ([2]), 
the rule space is condensed and a rule is now defined by 
the rule parameters (a, 7) and the threshold parameter 
n. 

As a first exemplary application, we consider scale-free 
graphs, generated by the incremental Barabasi-Albert 
model [17j, with a power law degree distribution, a prop- 
erty often found in real-life networks [181 [19] . Due to their 
pivotal topological property, the existence of hubs, scale- 
free graphs have been used frequently as model graphs. 
They have also been used as a starting point to investi- 
gate the relation of degree-degree correlations in complex 
networks [131 ISQl Elj • Here we want to study how degree 
correlations in a scale-free graph affect its ability to gen- 
erate complex patterns. We implement all resulting rules 
on randomized, hierarchized and anti-hierarchized vari- 
ants of scale-free graphs with 200 nodes and 400 links. 
Hierarchization and anti-hierarchization means the grad- 
ual randomization towards positive and negative degree- 
degree correlations, respectively [13]. For each graph 
type, we calculate the entropy signatures, given by the 
Shannon entropy and word entropy of the emerging pat- 
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FIG. 2: Number of complex patterns (solid colored lines) and 
number of rules (dashed line) for regular graphs of different 
networks size against degree d. The absolute number of possi- 
ble rules increases due to the enlarged possible set of n values. 
Strikingly, the number of single threshold rules with complex 
patterns is maximal at d = 4. 




0.5 0.0 0.5 

A Shannon entropy 

FIG. 3: Entropy signature difference plot for 25 scale- free 
graphs with 200 nodes and 400 links. Both positive and neg- 
ative entropy differences (A^ = Ey^ — E^) occur for the com- 
parison of hierarchized vs. randomized graphs. 
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terns, for all rules. Figure |3] shows the entropy signature 
difference plot, where {S^ W) of the randomized graphs 
(R) has been subtracted from the entropy signature of 
the hierarchized graphs (H), AE = E^l — £^r. Here, E 
stands for S and W respectively. Most rules are insen- 
sitive to degree-degree correlations. Positive and nega- 
tive entropy signature differences occur preferentially for 
{a,P) = (+,-) or = (-,+) with k e [0.3,0.6]. 

Notably, rule (+, — ) is a condensed form of the topology- 
sensitive ECA rule 108, appearing in Figure [l] For the 
anti-hierarchized graph with negative degree-degree cor- 
relations, a similar picture emerges (data not shown). 
As a second example, we consider the topology of 
metabolic networks, which abstracts the wiring architec- 
ture of the set of enzyme-catalyzed reactions in a spe- 
cific species. Substrate graphs, where nodes represent 
metabolites and links represent a reaction between con- 
nected substrates can be generated from genomic data 
[22] . In [12] , we recently studied the impact of the topol- 
ogy of metabolic networks on a specific dynamics. There 
we implemented and studied only a single rule, namely 
(+, — ) as a dynamic probe and interpreted the enhanced 
regularizing capacity of real networks compared to ran- 
domized null models as a possible topological contribu- 
tion to the reliable establishment of metabolic steady- 
states and to the effective dampening of fluctuations. To 
show that the results presented in [12 are valid over 
the whole range of dynamics discussed in the present 
paper, we now implement all possible rules of the form 
(a, 7) on substrate graphs and analyze the entropy sig- 
nature differences. Figure |4| shows the results for Homo 
sapiens^ Saccharomyces cerevisiae^ Escherichia coli^ and 
Bacillus subtilis. We find that while the entropy sig- 
natures of most rules do not discriminate between real 
and randomized topologies, a few rules are topologically 
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FIG. 4: Entropy signature difference plot for the metabolic 
networks of H. sapiens (N = 625, M = 779), S. cerevisiae 
(448, 564), E. coli (563, 709), and B. subtilis (501, 612). N 
denotes the number of nodes, M the number of links in the 
largest connected component of the substrate graphs of the re- 
spective species. Notably, predominantly negative differences 
AS and AW occur. 



sensitive. These rules comply with (a,/3) = (+, — ) or 
{a,p) = (-,+) while k G [0.35,0.55]. For ah these 
rules, the entropy signature of real graphs is significantly 
smaller compared to the null model topologies. This is 
also true for hierarchized and anti-hierarchized null mod- 
els, as well as for all other species investigated in [12]. We 
believe that the application of dynamic probes is a par- 
ticularly helpful tool for studying dynamical constraints 
imposed by topology. 

Discussion — Our formalism can be used to describe 
outer-totalistic CA and isotropic RBN rules in a com- 
mon framework. It allows the comprehensive discussion 
of previously introduced rule sets on diverse topologies, 
like the selection of Boolean rules presented in [7] or vari- 
ations of the majority rule as used in [6l[8]. It moreover 
formalizes previous attempts to generalize CA to graphs 
[26l[27], and is easily extensible, e.g. by introducing more 
than just one threshold parameter or by using a larger 
state space. With the presented framework, the often 
huge rule space of a discrete dynamical system can be in- 
tuitively parametrized and systematically analyzed. The 
finding of symmetric rules, for example, helps comple- 
menting specific CA classes. The set of rules exhibit- 
ing power law spectra, as introduced in [24], can thus 
be completed with the corresponding symmetric coun- 
terparts. Also, many coarse-graining transitions between 
CA rules, as presented in [28], can be immediately un- 
derstood with a symmetry rational. Specific rule gener- 
alizations, as discussed and analytically analyzed in [25 
may be reconsidered from the more general perspective 
provided in this letter. As a specific example, the ap- 
plication of ECA rule 22 to arbitrary graphs, stated as 
an open question in [25j , is straightforward with the pre- 
sented formalism. Limitations arise as soon as individual 
node characteristics are to be taken into account. Still, 
the isotropic subset of canalyzing Boolean rules, as dis- 
cussed in can be represented with our approach. 
Table [T] shows some examples of symmetric rules in our 



TABLE I: Examples of symmetric outer-totalistic CA as de- 
fined in Eq. |T]), corresponding ECA and RBN rule numbers, 
as defined in [2] and [i], respectively, and references where 
these rules have been discussed previously. 



Outer-totalistic CA 


ECA 


RBN 


Ref. 


(0,0,0)^(1,1,1) 


^ 255 


1 ^ 16 




(1,0,0)^(1,1,0) 


5 ^ 95 


2^8 


m 


(1,-,+) ^ (+,-,0) 


22 ^ 151 




[21 [25] 


(0,1,0)^(1,0,1) 


90 ^ 165 


7^ 10 


m 


(+,+,-) ^ (-,+,+) 


108 ^ 201 




[inidiKii] 


(+,!,-) ^ (-,0,+) 


126 ^ 129 




[21 [25] 


(0,0,1)^(0,1,1) 


160 ^ 250 


9 ^ 15 


[S] 


(0,1,+) ^(+,0,1) 


164 ^ 218 




m 


(0,+,l) 


232 
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formalism, the corresponding EGA and RBN rule num- 
ber and references where these rules have been previously 
applied. 

An analysis of topologically sensitive rules with analyti- 
cal tools as developed in [29] or the recently introduced 
basin entropy [30^ may reveal state space changes as- 
sociated with topological modifications. Such analyses 
can elucidate dynamic properties also relevant for regu- 
latory dynamics of biological networks, which have been 
successfully modeled with CA approaches [313 |32l [33] . 
From this perspective, our framework provides a means 
to comprehensively study the sensitivity of a system to 
topological perturbations and associated rule space mod- 
ifications. 
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